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Using ah initio density functional theory we explore the behavior of thin layers of metallic d* 
SrFeOs confined between the dielectric SrTiOs in a superlattice geometry. We find the presence 
of insulating SrTiOa spacer layers strongly affects the electronic properties of SrFeOs; For single 
SrFeOs layers constrained to their bulk cubic structure, the Fermi surface is two-dimensional, nested 
and resembles the hole-doped superconducting cuprates. A Jahn- Teller instability couples to an 
octahedral tilt mode, however, to remove this degeneracy resulting in insulating superlattices. 
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Progress in the layer-by-layer growth of transition metal 
oxide thin films [TJ [2] motivated the intriguing recent sug- 
gestion that oxide heterostructures engineered to have 
band structures close to those of the high- Tc cuprates 
could yield new superconductors [3]. Strained LaNiOs lay- 
ers separated by inert spacers such as LaAlOs or LaGaOa 
were proposed as a promising trial system, with strain 
lifting the degeneracy of the single Ni'^^ eg electron and 
layering providing quasi two-dimensionality. Indeed, first- 
principles electronic structure calculations 31 on strained 
Ni-based oxide superlattices have found that, with careful 
choice of strain and local- interface chemistry [5], 2D Fermi 
surfaces resembling those of the hole-doped cuprates ^ 
can be obtained. Superconductivity in such superlattices 
remains to be observed experimentally, however, possi- 
bly because competing instabilities such as charge- or 
orbital-ordering are enhanced by the reduced dimension- 
aUty. Although high-T^ superconductivity tends to occur 
in materials with large electronic fiuctuations and in close 
proximity to electronic, structural or magnetic phase tran- 
sitions O m m [TU] , the detailed role of such instabilities 
remains unclear [TTJ [12] . 

With these factors in mind, we examine how structural 
confinement and lattice instabilities modify the electronic 
structure of SrFeOs in SrTiOa/SrFeOa superlattices. Our 
motivation for this choice of system is three- fold: First, 
like Ni^+, Fe'*"'" has a single degenerate Cg electron in an oc- 
tahedral environment that dictates its low energy physics. 
The restriction on Fe-based compounds — traditionally 
dismissed from consideration for superconductivity be- 
cause of their robust magnetism — has been lifted due to 
the recent discovery of superconducting Fe pnictides [13' . 
Second, bulk SrFeOs is metallic with p-type conductivity 
(like the doped cuprates) and is proximal to multiple insta- 
bilities: It manifests a long-wavelength spin density wave, 
but neither Jahn- Teller distorts nor charge orders, even 
though both possibilities are suggested by its chemistry. 
Finally, unlike the structurally inert LaAlOa in the nicke- 
late superlattices, SrTiOs is a highly polarizable dielectric 
which can couple to electronic or structural distortions 
[Ti] in the SrFeOa layer. 



Using first-principles density functional theory (DFT) 
within the local-spin density approximation (LSDA) plus 
Hubbard U method, we calculate the structure and elec- 
tronic properties of (SrTi03)„/(SrFe03)m superlattices. 
We focus on (?) the evolution of the 2D band structure 
with thickness of the dielectric and metal (n, to = 1 or 
3), and (iz) how competing structural and electronic in- 
stabilities manifest in the superlattices. We find that 
2D confinement from the superlattice periodicity along 
the growth (z) direction yields low-energy physics that 
are primarily derived from hybridized Fe dx^-y^ and O 
2p orbitals. The idealized high-symmetry superlattice 
structure has a strongly nested 2D Fermi surface that is 
similar to that of the parent superconductor La2Cu04 [TS]. 
Strong electron-lattice instabilities, enhanced by 2D con- 
finement, transform the metastable metallic structure into 
a lower symmetry insulating superlattice. Understanding 
the competing nature of these often neglected lattice dis- 
tortions is crucial for engineering electronic structures of 
oxide superlattices and interfaces. 

Our LSDA-|-[/ DFT calculations are performed using 
the Vienna ah initio Simulation Package (VASP) |TF1 ITT] . 
We follow the Dudarev approach [18] and include an effec- 
tive Hubbard term U^s = C/ — J of 6 eV to treat the Fe Sd 
orbitals [3D]. This method and value of U^e gave good re- 
sults in earlier first-principles calculations for bulk SrFeOs 
and related iron oxide compounds [THl HO]- We construct 
the superlattices by stacking five-atom perovskite units 
along the z-direction [see Fig. [ijd)], and constrain the in- 
plane lattice parameter (xy-plane) to that of cubic LDA 
SrTiOa (a = 3.86 A) to simulate growth on a SrTiOa 
substrate. While the lattice mismatch between SrFeOa 
and SrTiOa is small (the theoretical mismatch is 1.8%, 
with ferromagnetic LSDA-|-t/ SrFeOa having a = 3.79 A), 
we show below that it is sufficient to impose an epitaxial 
crystal field (Aecf) that partially lifts the Fe Cg orbital 
degeneracy. We then relax the length of the c-axis and 
the internal degrees of freedom along the z-direction. The 
in-plane periodicity imposed by this choice of reference 
structure does not permit cell-doubhng structural distor- 
tions such as octahedral rotations or orbital ordering: This 
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FIG. 1: (Color) The DOS for (a) cubic SrTiOs (a ^ 
(b) cubic FM SrFeOg (a = 3.79 A), and (c) the high symmetry 
PA/mmm (SrTi03)3/(SrFe03)i superlattice shown in (d). 



highly-symmetric superlattice (space group PA/mmm) 
resembles the structural constraints often considered to 
be imposed in heteroepitaxial thin films. We relax this re- 
striction later to study electron-driven lattice instabilities. 
Ferromagnetic (FM) order is imposed on the Fe sites in 
all calculations, and is theoretically [20] found to be the 
lowest energy coUinear ordering for bulk SrFeOa. 

We begin by studying the (SrTi03)3/(SrFe03)i super- 
lattice and plot in Figure [T] our calculated local densities- 
of-states (DOS) for (a) cubic SrTiOg, (b) cubic SrFeOg, 
and (c) the superlattice. Our results for the bulk com- 
pounds are consistent with the literature: fn SrTiOs we 
obtain a ^^2 eV band gap between an O 2p valence band 
and a Ti 3d conduction band. SrFeOs is nearly half- 
metallic with a calculated magnetic moment of 3.8 /is 
per Fe atom, consistent with a high-spin electronic 
configuration. The Fermi level (dashed line at eV) lies 
in a region of majority spin Fe eJ-O 2p^ hybridized or- 
bitals and in the O 2p^ valence band. The electronic 
structure of the superlattice is close to a superposition 
of its constituents, except for a reduction in the contri- 
bution of the O 2p states at Ep making the superlattice 
half-metallic. While the use of the Hubbard U method in 
DFT calculations is known to drive bands toward integer 
filling, the half-metallicity is robust for [/ = to 8 eV. 

We plot the band structure of the PA/mmm superlat- 
tice in Figure [2]^a,b) using the "fat-band" method [21 in 
which the magnitude of the projection of each Bloch state 
onto a particular set of atomic orbitals is represented by 
its line width: Panel (a) shows the equatorial O 2px,y-Fe 
3dj;2-y2 states and (b) shows the O 2p^-Fe 3dz2 orbitals. 
Consistent with the DOS, we find partial occupation of 
the majority dx2_y2 and d^^ Fe states at Ep with the 
energy of the dxi_yi at F 0.50 eV lower than that of 
the dz2 orbital. The energy difference at F is explained 
from compression of the apical (1.88 A) over the in-plane 
(1.93 A) Fe-O bond length from the imposed Aecf 
Interestingly, the apical Fe-O distance is shorter than the 
value (1.90 A) that we obtain when we force hulk SrFeOa 
to have the SrTiOs in-plane lattice constant, suggestive of 



an additional spontaneous Jahn- Teller-like distortion (we 
examine this later) . Despite the substrate elongation of 
the equatorial Fe-O bonds, the bands of dx2-y2 character 
do not fully split from the bands to produce a single 
band crossing Ep] in fact, both the top and the bottom 
of the Cg band are dominated by dx2^y2 states. The 
dz2 band can likely be further destabilized to achieve an 
identical cuprate-like electronic structure with additional 
strain [3] and/or chemical modifications [J. 

We now investigate the character of the 2D bands and 
examine how the electronic structure aX Ep responds to 
changes in the dielectric thickness by comparing the full 
band structures of the (SrTi03)3/(SrFe03)i [Fig.|2](a,b)] 
and (SrTi03)i/(SrFe03)i [(c,d)] superlattices. As ex- 
pected, the bandwidth of the dx2-y2 derived states along 
F-A is nearly the same in the two superlattices due to 
identical in-plane structural parameters. The dispersion 
of the ^2-2 bands is markedly different, however, with the 
less confined n — l,m — 1 superlattice having a larger 
bandwidth than the more confined n — 3,m — 1 along 
the M-A lines. Increasing dielectric thickness {n — 1 ... 4) 
shows that the dispersive dx2__y2 bands saturate at ~2 eV 
in width while the d^2 states are almost dispersionless, 
although never fully split. In contrast, increasing the 
thickness of the SrFe03 layers in the superlattice pro- 
duces several partially occupied eg bands crossing Ep, 
and the band structure (not shown) resembles a regular 
metal. We conclude that a cuprate-like band structure is 
unlikely in superlattices containing multiple ferrate layers. 

The 2D confinement is discernible in the Fermi surface 
of the n = 3,m = 1 superlattice shown in Fig. [2je,f). 
The additional band crossing Ep in the kz ~ plane. 
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FIG. 2: (Color) Electronic band structure plots for the n = 
3, m = 1 [panels (a) and (b)], and n = 1, m = 1 [panels (c) and 
(d)] superlattices. Panels (a) and (c) show the majority spin 
"fat-bands" (blue) derived from O 2px,y-Fe 3d^2_y2 states, 
whereas (b) and (d) show the O 2pz-Fe 3dz2 states. The 
Fermi surface at fcz = (e) and in the full (f) BZ shown in 
(g) consists of two bands connected by the Q — (1,1,0)^ 
interband nesting vector (arrow). 



3 



not seen in the n = 1, m = 1 (c,d) super lattice nor bulk 
SrFeOs lattice matched to SrTiOa, produces squared- 
cylindrical arrays around the M-point consistent with the 
four-fold symmetry of the lattice. An unusual inward 
bowing along the F-X direction occurs with corrugations 
along kz similar to over-doped cuprates. The band cur- 
vature indicates electron and hole sheets at T and M, 
respectively, similar in character to the superconducting 
cuprates \T5[ [^ . Interestingly, the two Fermi sheets 
nearly intersect at half the F-M distance producing a 
nesting vector Q = (1, 1,0)^ connecting them. Such a 
vector can result in a charge or spin density wave, or for 
strong coupling of the electronic system to the lattice, 
a symmetry lowering structural distortion. These two 
sheets are largely insensitive to the calculation details: 
changing U weakly affects the band crossings at Ep since 
the electronic structure is mainly modified by the con- 
fining SrTiOs layers. Spin-orbit interactions also do not 
alter the FS nesting. 

We now investigate possible structural instabilities in 
the (SrTi03)3/(SrFe03)i superlattice, and their influence 
on the FS degeneracy. We identify likely M-point insta- 
bilities (consistent with the nesting vector) by choosing 
structural distortions that connect the high symmetry ref- 
erence phase to low symmetry structures through atomic 
distortions that maintain a direct group-subgroup rela- 
tion. We choose these distortions to be irreducible 
representations (irreps) of space group PA/mmm and 
consider only M^, for i — 1...5. (We ignore M~ ir- 
reps with antisymmetric distortions under inversion, and 
follow the notation of Miller & Love throughout |24j.) 
Fully relaxed ab initio structures are then reduced into 
combinations of these irreps by performing a symmetry 
mode analysis pS] that makes accessible the local atomic 
displacements. We keep FM spin order fixed as above to 
isolate the electron-lattice coupling. 

The active irreps we examine (Fig. [s]) affect either the 
equatorial Fe-O bond lengths (crystal field) or 0-Fe-O 
bond angles {dpa bandwidth). Irreps and are pla- 
nar breathing and stretch modes of the Fe04 plaquettes, 
respectively, and are anticipated to strongly modulate the 
eJ-O 2p hybridization at Ep- The M^ mode produces 
two unique Fe04 plaquettes arranged in a 2D checker- 
board manner and gives rise to charge disproportionation 
in bulk isoelectronic CaFeOa In contrast the ir- 
rep, or Jahn- Teller mode, creates two long and two short 
equatorial Fe-O bonds and favors orbital polarization. 
The remaining irreps are collective distortions of the oc- 
tahedral units that produce the common tilt and rotation 
patterns found in perovskite oxides. These modes affect 
the dpa bandwidth through deviations in the Fe-O-Fe 
bond angle away from ideal 180°: M;^ consists of rota- 
tions of the in-plane oxygen atoms about the z-asds; 
produces a bending distortion of the 0-Fe-O bond; and 
the 2D irrep M^ produces two possible tilt patterns of 
the planar oxygen atoms: M^(a,0) giving a pattern with 
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FIG. 3: (Color) Electronic responses from SrFeOs (orange), 
SrTiOa (blue), and the superlattice (green) to irreps 
of space group PA/mmm: M~[ breathing mode, stretch 
mode, Mg" rotation mode of the FeOe octahedra around the 
z-direction, bending of the Fe-O bonds (purple), and two 
different tilt patterns M^(a, 0) and M^(a,a) in the xy-plwas. 



space group symmetry Pmna or M^(a,a) with Cmma. 

We begin by "freezing" each irrep as a function of mode 
amplitude into the equatorial oxygen atoms coordinating 
Fe (see Supp. Materials). At the same time, we keep 
the remaining atoms in the superlattice clamped to the 
reference PA/mmm positions to isolate the response of 
the ferrate layer — displacements of other atoms in the 
superlattice are obtained later through full structural 
optimization. The energy change due to the optimal 
mode amplitude for each irrep frozen in the ferrate layer 
is shown in Fig. [3] with respect to the high symmetry 
PA/mmm superlattice (energies per five-atom perovskite 
cell). Irreps M^, M^ and Mj^ lower the energy of the 
superlattice. In contrast they increase the energy of hulk 
SrFeOs lattice matched to SrTiOa, indicating that the 
superlattice geometry makes the eg orbital degeneracy 
more susceptible to structural distortions. The remaining 
irreps which require collective cooperation of all oxygen 
octahedra increase the superlattice energy when the di- 
electric is clamped to the reference configuration. Next, 
we fully relax the SrTiOa with the ferrate atoms fixed to 
the distortions above to isolate the effect of the dielectric 
response. Now every irrep except MJ is found to be en- 
ergy lowering: the large energy gain now observed for the 
M^ modes suggests that a 3D tilt pattern is preferred — 
this result is likely driven by the fact that most of the 
superlattice is SrTiOa which itself has a large M^-like 
mode instability. 

We now relax both constituents under the symmetry 
constraints imposed by the irreps to find the lowest en- 
ergy atomic configurations. The largest energy-lowering 
distortions are found to be the Jahn- Teller and tilt modes. 
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FIG. 4: Band structure for (a) the Pi/mmm [n — 3,m = 1) 
superlattice as in Fig. [2|a) with the lattice parameter doubled 
in-plane to allow comparison with panels (b)-(d) which include 
additional structural distortions: (b) M]^, (c) MJ, and (d) 
M^(a,0); and the ground state P2\/c structure (e). 

Both irreps and M2', open energy gaps at the Fermi 
level in the superlattice [Fig. |4]jb,c)], while identical dis- 
tortions in bulk SrFeOs maintain its metallicity. The 
tilt modes M^(a,0) and M^(a,a) are more unstable than 
the rotation irrep with M^(a,0) almost as energet- 
ically favorable as the Jahn- Teller distortion. Yet, due 
to the weaker influence of the bond angle on the dpa 
bandwidth, none of the bond angle distortions gap the 
FS. There is a 0.24 eV bandwidth reduction, however, for 
irrep shown in Fig. |4|d) . Similarly, the mode also 
gaps the Fermi surface — although not an energy lowering 
distortion — in doped La2Cu04, whereas the octahedra 
tilt and rotation modes, which are found in the cuprates, 
show minor effects [21 . 

Finally, we calculate the fully optimized structure start- 
ing from a combination of the and irreps. Our 
ground state structure is 33.5 meV lower in energy (per 5- 
atom unit cell) than the reference superlattice. It contains 
Fe04 plaquettes with bond lengths of 1.86 and 2.01 A 
due to the Jahn-TelUer mode thats makes the superlattice 
insulating [0.58 eV gap. Fig. |4|e)]. The octahedral insta- 
bilities largely condense in the SrTiOs layers with reduced 
interaction in the ferrate plane: The in-plane Fe-O-Fe 
bond angle is reduced by 6.89° while a 10° reduction 
occurs across the Fe-O-Ti angle. By symmetry decom- 
posing the structure as a combination of the PA/mmm 
irreps, 0.05ir+ + 0.048r+ + 0.21M+ + 0.13M+ + 0.97M+, 
we find the primary order parameters and drive 
the P4/mmm (a°,a°,c°) P2i/c {a+,b~,b~) metal- 
insulator transition, with the associated change in octa- 
hedral tilt patterns. Since none of these distortions are 
favored in bulk SrFeOa, we attribute the electron-lattice 
coupling enhancement to the confinement effects imposed 
by the superlattice geometry. 

We have demonstrated that in a superlattice geome- 
try, confinement from dielectric spacer layers combines 
with a strain-induced epitaxial crystal field (Aecf) to 
modify the energies and dispersions of apical d^^ and 
in-plane Fe dx2_y2 orbitals. When we simulate the high 
symmetry cubic structures of room temperature SrTiOs 
and SrFeOa partial planar electron localization yields 
a 2D Fermi surface that resembles the superconducting 
cuprates. Low-energy M-point instabilities compete with 
nesting on the Fermi surface to make the superlattices 



proximal to multiple competing phases. The lattice insta- 
bilities are enhanced by the dielectric SrTiOa layers, there- 
fore inert spacer layers without polarizable ions might be 
more favorable for superconductivity. We hope that our 
finding of competing structural ground states with distinct 
electronic properties motivates experimental investigation 
of SrFeOa/SrTiOs heterostructures using external electric 
and magnetic fields: Such probes could tune between the 
itinerant and localized electronic states in analogy with 
the parent Mott insulating cuprates. 
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